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Abstract 

As a first result of two new partial- wave analyses, one of the pp and another one of the np scatter- 
ing data below 500 MeV, we report a study of the long-range chiral two-pion exchange interaction 
which contains the chiral coupling constants ci, C3, and C4. By using as input a theoretical value 
for ci we are able to determine in pp as well as in np scattering accurate values for C3 and C4. 
The values determined from the pp data and independently from the np data are in very good 
agreement, indicating the correctness of the chiral two-pion exchange interaction. The weighted 
averages are C3 = — 4.78(10)/GeV and C4 = 3.96(22)/GeV, where the errors are statistical. The 
value of C3 is best determined from the pp data and that of C4 from the np data. 
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Introduction 

It is beyond doubt that the longest-range strong two-nucleon (AW) interaction is the 
one-pion exchange (OPE) force. Despite more than 50 years of research, the nature of 
the medium-range AW interaction is not so well understood. What seems clear is that 
it contains: (i) A strongly attractive central force, (ii) an isospin-dependent tensor force 
opposite in sign to OPE, and [in) a rather strongly attractive spin-orbit force. It was 
discovered in the early sixties that all these features follow naturally from the exchange of 
scalar and vector mesons, which led to the development of the one-boson exchange (OBE) 
model of the WW interaction. The role of the two-pion exchange (TPE) interaction and its 
interplay with the exchange of heavy mesons that decay into two pions has for a long time 
remained elusive. 

In recent years, however, the situation has improved. The derivation of at least the long- 
and medium-range nuclear forces can be formulated in a model-inclCT)endent manner by a 

systematic expansion of the chiral Lagrangian of QCD Biiiiiaia. 

In particular, 

the long-range TPE interaction can be derived unambiguously, where the effects of the 
exchange of broad heavy mesons are incorporated in effective low-energy chiral coupling 
constants. Most importantly, chiral symmetry and its breaking are correctly implemented 
in this approach. 

In Ref. we studied this long-range chiral two-pion exchange (xTPE) interaction in 
an energy-dependent partial-wave analysis (PWA) of the proton-proton (pp) scattering data 
below 350 MeV. The presence of xTPE in the long-range pp force was demonstrated, and 
the chiral coupling constants C3 and C4 were determined from the pp data. In this paper 
we address the question whether the same xTPE force allows also a good description of 
the neutron-proton [np) scattering data below 500 MeV. Moreover, we present new, precise 
determinations of the chiral coupling constants C3 and C4 from the pp and np data separately. 
Accurate values of these chiral coupling constants are an important input in calculations of, 
for instance, the two-pion exchange three-nucleon force 0]. 



Partial-wave analysis 

Because of the high quality of the pp data base all the pp phase shifts with orbital angular 
momentum i < 4 can be determined accurately in an energy- dependent PWA of the pp data 
below 350 MeV. An analysis of the np data, however, is much more difficult, because not 
only the 1 = 1 phase shifts but also the 1 = phase shifts contribute. Moreover, the np 
data base, while extensive, is by far not as accurate and varied as the pp data base jij. In 
a PWA of only the np data below 350 MeV it has always been impossible to determine all 
the important np phase shifts. Therefore, the standard practice has been to take the 1 = 1 
phase shifts, with the exception of the ^5*0 phase shift, from the pp PWA, with or without 
corrections for the Coulomb interaction and/or the tt^-tt^ mass difference in OPE. Since it 
has long been known that there is a sizable charge-independence breaking (GIB) in the ^Sq 
phase shifts, the ^Sq np phase shift is always fitted independently of the ^5*0 pp phase shift. 

This approach to np PWA was also followed in the past by the Nijmegen group. In 1993 
the results of the first Nijmegen pp and np PWA's below 350 MeV were published in Ref. jlol |. 
An attempt at that time to extract all the important np phase shifts, both 1 = and 1 = 1, 
from the np data base below 350 MeV failed, although it was possible to determine the 
np phase shifts when the 7 = 1 waves for i > 1 were taken over from the pp PWA93 
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It has been customary to perform AW PWA's without inelasticities up to 350 MeV, 
ahhough pion production starts aheady at 280 MeV. It can be shown that the inclusion 
of inelasticities in the pp PWA below 350 MeV improves the Xmm slightly. Already some 
time ago [1^ the Nijmegen pp PWA was extended to energies far above the pion-production 
thresholds, with the inclusion of inelasticities. When, in 1994, the np PWA was extended 
to 500 MeV, it turned out to be possible, for the first time, to determine uniquely all the 
important np phase shifts, both 7 = and 7 = 1, from the np data alone [13[. Such a 
separate PWA of the np data, without input from the pp PWA for the 7 = 1 waves, is 
in principle more model independent. A comparison between the phase shifts from such 
an independent np PWA and the corresponding phase shifts from the pp PWA provides 
information about possible CIB in the 7 = 1 waves. 

The Nijmegen energy-dependent PWA's can be used as a tool to study the long-range 
NN interaction |7]. The long-range forces are included exactly, in order to ensure that the 
partial-wave amplitudes acquire the proper fast energy dependence from the nearby left-hand 
singularies due to these long-range forces, while the short-range interactions (more remote 
left-hand singularities), responsible for a much slower energy dependence, are parametrized. 
This strategy is implemented by solving the Schrodinger equation with an energy- dependent 
boundary condition (BC) at some r = b and for r > b the long-range NN interaction. This 
long-range force contains the electromagnetic interaction (i.e., in the pp case the improved 
Coulomb the magnetic- moment [l^, and the vacuum-polarization fl6^ interactions, 
and in the np case the magnetic- moment interaction Il5l'l . the OPE interaction, in the np 
case also the pion-photon (n-'j) exchange interaction [17|, and the long-range part of the 
xTPE interaction [7|. The BC is parametrized as an analytic function of energy, and the 
parameters, representing "short-range physics," are determined from a fit to the data. The 
option also exists to fit simultaneously some of the parameters in the long-range interactions, 
viz. the pion-nucleon coupling constants 3, IS 2^ and/or the chiral coupling constants q 



(z = 1,3,4) in xTPE [7]. 

The new pp and np PWA's that we discuss here will be referred to as xPWA03. They 
differ from the old PWA93 in several aspects: (i) The energy range is extended from 350 
to 500 MeV. Instead of the 1787 pp data and 2514 np data in PWA93 we now have 5109 
pp data and 4786 np data 21]. (ii) All the np phase shifts can be determined from the np 
data alone, instead of taking the 7=1 phases from the pp PWA and correcting them. [Hi) 
Inelasticities are taken into account, (iv) For r > 6 a different non-OPE strong interaction is 
taken. In PWA93 the heavy-meson exchanges of the Nijmegen soft-core OBE potential 
were used. Motivated by the success of Ref . , where an excellent description of the high- 
quality pp data base below 350 MeV was obtained, we use in xPWA03 the xTPE potential. 
(v) A minor difference between PWA93 and xPWA03 is that we take here b = 1.6 fm, while 
in PWA93 b = 1.4 fm was used; in Ref. [tI we used both b = 1.4 fm and b = 1.8 fm. 

Details of xPWA03 (data, phase shifts, etc.) will be presented elsewhere ji^, here we 
focus on testing the long-range xTPE interaction in the pp and np systems below 500 MeV. 



Chiral two-pion exchange potential 

The xTPE potential can be derived by a systematic expansion of the effective chiral 
Lagrangian [11,13, IE B S S Q- The form that is appropriate for use in the relativis- 
tic Schrodinger equation and that is consistent with our choice of including the minimal- 
relativity factor M/E in the OPE potential is specified in Ref. (see also Ref. ^). 
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The leading-order xTPE potential consists of the static planar- and crossed-box TPE 
diagrams, calculated with the derivative (pseudovector) NNn Lagrangian, and the triangle 
and football diagrams with the nonlinear NNnn Weinberg- Tomozawa (WT) seagull vertices. 
It contains isospin-independent spin-spin and tensor terms and an isospin-dependent central 
term. 

In subleading order, next to nonstatic corrections to the planar and box diagrams, ad- 
ditional triangle diagrams appear which contain three new NNtttt interactions The 
corresponding chiral coupling constants are denoted by q {i = 1,3,4). (Unfortunately, they 
are not scaled to obtain dimensionless numbers and their values are conventionally given in 
GeV~^.) They are defined by the following terms in the chiral Lagrangian density: 



£ = -N 



SciD-'mlTryF^ + AcsD^-D'' 



+ 2c4a^.r-5^xD"] iV , (1) 

where F^^ ~ 185 MeV is the pion decay constant, D = 1 -\- tt^/F^, and the chiral- covariant 
derivative of the pion field vf is = D^^d^Ti/ F.^. The C3- and C4-terms are manifestly chiral 
invariant. The Ci-term, which is proportional to m^, violates chiral symmetry explicitly 
and is related to the much-discussed pion-nucleon sigma term ^|. Using the rationalized 
pseudovector pion-nucleon coupling constant p the relation reads [25^ 



ci = - 



Ami 16 



s 



(2) 



where m.,^ = 138.04 MeV is the average pion mass, and = is the scaling mass 
conventionally introduced to make / dimensionless. Eq. (E)) holds in order 0{q^) in the 
chiral expansion in small momenta q and the pion mass ^]. An additional C2-term is not 
given in Eq. (jT)), since it does not contribute to the xTPE potential to subleading order. 
However, it does contribute to the isoscalar ttN scattering amplitude at the same order as 
ci and C3. 

In subleading order the xTPE potential gets contributions to the central, spin-spin, ten- 
sor, and spin-orbit potentials (c/. Table 1 in Ref. Q). Important components are: (i) A 
strong isospin-independent central attraction due to the Ca-term, (u) an isospin-dependent 
tensor force opposite in sign to OPE due to the C4-term, and (iii) an attractive isospin- 
independent spin-orbit force from nonstatic terms of the planar- and crossed-box diagrams. 
The values of the Cj's are not fixed by chiral symmetry and must be determined from the 
experimental vrA^ or NN scattering data. 

The long-range xTPE potential derived in the framework of the effective chiral La- 
grangian is completely model independent. Any dynamical model [IH for the TPE 
NN interaction, containing e.g. the e (or "a") and g mesons, the pomeron, and/or A^- and 
A-isobars, has to reduce to this form for large r. These models should also predict values 
for the Cj's consistent with the determinations from the vrA^ and AW scattering data. 

The breaking of charge-independence due to the tt'^-tt" mass difference in the OPE poten- 
tial is taken into account exactly, as it was already in PWA93 H^- In thevTPE potential 
we include the terms linear in the 7r"'"-7r° mass difference, following Ref. 2^|. One charge- 



independent pion-nucleon coupling constant i20| is used in both the OPE and the xTPE 
potentials. In the long-range interaction for r > b only the chiral coupling constants q 
{i = 1,3,4) remain then to be determined. 
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Results 



In our previous study [7] of xTPE it turned out that ci could not be determined accurately 
from the pp data base below 350 MeV. When we fitted Ci, C3, and C4 simultaneously, we found 
Ci = — 4.4(3. 4)/GeV, where the error is statistical. A strong correlation was obtained be- 
tween the values of Ci and C3. Therefore, we used Eq. (0) to fix the value of Ci. Assuming that 
the Sigma term has the low value a = 35(5) MeV ^2^,12^], Eq. © gives Ci = -0.76(7)/GeV, 
where the error is theoretical. We used the central value ci = — 0.76/GeV as input in the 
PWA, and determined in Ref. the values C3 = -5.08(28)/GeV and C4 = 4.70(70)/GeV, 
where the errors are statistical. One notes that from the pp data below 350 MeV the value 
of C3 could be extracted rather precisely, while C4 was pinned down less accurately. 

The value extracted for ci from the data below 500 MeV would also not be accurate 
enough to shed light on the value of the sigma term, since the statistical error for ci obtained 
in Ref. would have to be reduced at least by a factor of about 20. We therefore decided 
to use also here the value Ci = — 0.76/GeV as input value, and to determine C3 and C4 from 
direct fits to the pp data and independently also from fits to the np data. 

We analyzed 5109 pp data below 500 MeV using 33 BC parameters, and we reached 
Xmin = 5184.3. The optimal values for C3 and C4 and their (1 s.d.) errors as determined 
from these pp data are: 



where also the dependence on the NNtt coupling constant is displayed. The correlation 
parameter is g = —0.47. These values for C3 and C4 are consistent with and more accurate 
than those found in Ref. from the pp data below 350 MeV. The errors are statistically 
only. Systematic errors are difficult to assess and require further study. 

For np scattering we analyzed 4786 data below 500 MeV. In this case we needed 40 
BC parameters and reached Xmin — 4806.2. The chiral coupling constants, their statistical 
errors, and their dependence on the NNn coupling constant are in this np case: 



The correlation parameter is g = 0.22. In Fig. we show the results for C3 and C4, for 
P = 0.0755. Plotted are the positions of the x^-minima and the = Xmin + 1 ellipses in 
the (c3, C4) plane, both for the pp and the np case. (These ellipses, of course, are determined 
with optimalization of all the BC parameters.) 

The values of C3 and C4 determined from the pp and from the np data are in good 
agreement. The value for C3 determined from the pp data is more than twice as accurate 
as the value from the np data, while for C4 the situation is reversed: the value from the np 
data is twice as accurate as the value from the pp data. 

We also determined the weighted averages with errors. We get 



-4.78(11) + 80(/2 - 0.0755)J /GeV , 
3.92(52) + 260{f - 0.0755)1 /GeV , 



(3) 



C3 = [-4.77(22) + 100(/2 - 0.0755)J /GeV , 
C4 = [3.97(24) + 40(/2 - 0.0755)1 /GeV . 



(4) 



C3 = [-4.78(10) + 84(/2 - 0.0755)J /GeV 
C4 = [3.96(22) + 79{f - 0.0755)1 /GeV . 



(5) 
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FIG. 1: Ellipses of constant in the (03, C4) plane. Shown are the = Xmin + -'- ellipses in the 
pp PWA and in the np PWA. The centers of the ellipses correspond to the minima in x^- 

These are our best values, following from all the pp and np scattering data below 500 MeV, 
which amounts to a total of almost ten thousand NN data. In Table E] we list our results for 
C3 and C4, for p = 0.0755. 



Discussion and summary 

We have determined accurate values for the chiral coupling constants C3 and C4 from the 
pp and the np scattering data below 500 MeV. The values for the Cj's {i = 1, 2, 3, 4) can also 
be determined from PWA's of the irN scattering data by fitting the amplitudes predicted by 
(heavy-baryon) chiral perturbation theory (xPT) to the nN scattering amplitudes obtained 
from these PWA's. In the several such determinations (for a discussion of the status see 
Ref. the value of C3 is found to lie in the range between —4.70 and —6.19, and the 
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TABLE I: Comparison of the chiral coupling constants q {i = 1,3,4) (in units of 1/GeV) from 
different analyses. The vrA^ results correspond to analyses to order in xPT. The input values 

of the sigma term are in MeV. (For a discussion of the meaning of the errors, see the text.) 



Ref. 


a 


Cl 


C3 


C4 


ttN 


[31] 


45(8) 


-0.91(9) 


-5.16(25) 


3.63(10) 


ttN 


[32] 


40(8) 


-0.81(12) 


-4.70(1.16) 


3.40(04) 


pp 


m 


35(5) 


-0.76(7) 


-5.08(28) 


4.70(70) 


pp 


This work 


35(5) 


-0.76(7) 


-4.78(11) 


3.92(52) 


np 


This work 


35(5) 


-0.76(7) 


-4.77(22) 


3.97(24) 


NN 


This work 


35(5) 


-0.76(7) 


-4.78(10) 


3.96(22) 



value of C4 in the range between 3.25 and 4.12. The values that are found for cs and C4 
depend on the order in xPT of the amplitudes, as well as on the specific vrA^ PWA that 
is used. They also depend on what value is used for the sigma term, because this value 
fixes the value of Ci. In xPT the isovector 71 N amplitudes can be predicted more accurately 
than the isoscalar amplitudes, because in leading order the latter are zero. Therefore, C4 can 
probably be pinned down better than C3. The value of C3 is moreover strongly correlated 
with the values of Ci and C2. In Table Ewe also listed the values obtained for the Cj's in 
two vrA^ analyses (to order O(q^) in xPT) that assume, like us, an acceptably low value for 
the sigma term. (From Ref. j3l| we list only the q's corresponding to one of the fits with 
/2 = 0.076.) 

A problem with these determinations of the q's from the vrA^ scattering data is that they 
are not determined directly from the data, but from fitting to the amplitudes of existing 
PWA's that have no reliable errors. For instance, in several analyses the amplitudes of the 
about 25-year-old Karlsruhe-Helsinki dispersion analysis were used. The amplitudes of that 
PWA have no associated errors and, what is worse, are in disagreement with the modern- 
day ttN data base. That analysis produced the high value = 0.079 for the pion-nucleon 
coupling constant |2^. The resulting errors on the Cj's determined from the vrA^ data, 
therefore, do not reflect the statistics of the data base, but are essentially rather arbitrary 
estimates. 

Our results correspond to a long-range interaction that includes the leading and sublead- 
ing xTPE diagrams. Higher-order corrections for xTPE and the leading three-pion exchange 
diagrams have been calculated by Kaiser [s^ , and can in principle be included as well. Work 
along these lines is continuing. 

In summary, the long-range part of the xTPE potential was included in energy-dependent 
PWA's of the pp and the np scattering data below 500 MeV. Good fits to the data were 
obtained. In the np PWA all the phase-shift parameters could be determined without input 
from the pp system. We conclude that OPE plus xTPE provides a high-quality long-range 
strong two-nucleon interaction. Accurate values for the chiral coupling constants C3 and C4 
of chiral perturbation theory were obtained from the pp and the np data separately. The 
values agree very well with each other, and they are also in good agreement with the range 
of values obtained from pion-nucleon scattering amplitudes. We consider this agreement to 
be experimental evidence that the xTPE interaction, as predicted by chiral perturbation 
theory, is correct. 
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